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Abstract 

Using the General Cascade Program(GCP), the production and absorption 
of J/ij: in p-A and A-A collisions have been studied. Nucleon absorption 
mechanism and comover absorption mechanism are considered to investi- 
gate the J/ij: suppression. The results agree well with experimental data of 
J/ijj production, except for the data in Pb-Pb collision. 
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I. Introduction 

Suppression of J/ip production in high energy heavy ion colhsions was proposed as a 
promising signature for the formation of quark gluon plasma(QGP) by Matsui and Satz 
ten years ago This suppression effect was observed by NA38 collaboration later |@]. 
However, it has been found that J/ ip suppression exists also in p-A collisions where QGP 
formation is not possible |^. Since then, the source of the observed J/ip suppression has 
remained controversial. Recently the NA50 collaboration reported the anomalous J /if) 
suppression in Pb-Pb collision [^,0. Now the debate has focussed on whether or not the 
data can be consistently described by hadronic mechanisms 0-10]. So it is important to 
study first the J / suppression caused by hadronic mechanisms. 

There are several sources of J/ip suppression in hadronic matter. Such as the absorp- 
tion of J/ip in nuclear matter, the interaction of J/ip with the produced mesons (called 
comover), gluon shadowing in nuclei, energy degradation of produced J/ip. In this paper, 
we will discuss the effects of the first two hadronic mechanisms on J/ip suppression. 

General Cascade Program(GCP) is designed by Y. Pang. It is a Monte Carlo 
simulation code based on cascade method. Cascade method has been used extensively in 
studying relativistic heavy ion collisions. It is perhaps the only quantitative tool currently 
available which is capable to provide both the overall features and the specific properties 
of nucleus-nucleus collisions at extremely high energies. 

In principle, the cascade algorithm itself is rather simple. One could visualize most 
cascades as collisions of classical billiard balls traveling at relativistic speed. Most cascade 
codes, on the other hand, are often quite complex. This is because each code must 
include a large number of physical processes in order to generate a realistic nucleus- 
nucleus collision event. Additional assumptions is often required in order to simulate an 
actual nucleus-nucleus collision. There are also parameters introduced during the cascade. 
The final result of a cascade often depends on these assumptions and parameters. The 
physical processes in these models can be quite different. For most cascade based models, 
it is often difficult to extract the key physics ingredient without having to go over the 
entire source code. It is not easy to compare various models, even when good agreements 
between the model predictions and the experimental data in the observed spectra have 
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been reached. 

GCP is designed such that the cascade algorithm is separate from the complexities of 
physics models. So the problem of the cascade can be isolated and studied separately, and 
different physics models can be compared in a common framework. GCP itself is only a 
cascade model, a tool for building relativistic cascade models including various physical 
contents. This is the key difference between GCP and most other cascade codes. 

By now there are many discussion on J/ip suppression. However, only a few works 
started to treat the Monte Carlo simulation of J/ip production and absorption directly 



T^,[T3|. While [jT2[ claimed that the J /ip suppressions in p-A and A-A collisions, including 



Pb-Pb data, could be explained based on nucleon and comover absorption, ref [|T3|| showed 



that it is impossible to fit the abnormal J/ip suppression in Pb-Pb collision only based 
on hadronic absorptions. Besides, the influence of the lifetime of the wounded nucleon is 
not investigated in details before. This is the reason why we use the code GCP to study 
J/ip suppression. We intend to study the influence of the lifetime of wounded nucleon 
on the strength of comover absorption. Our results show that it is not possible to fit the 
Pb-Pb data within reasonable choices of parameters. Our work is divided into two parts. 
First, the longitudinally excited wounded nucleon model(LEWNM) is added to GCP. We 
use it to handle the nucleus-nucleus interaction. Next, the process of J jip production and 
absorption in p-A and A-A collisions in SPS energy is simulated using GCP. The results 
agree well with most experiment data of J / ip production in p-A and A-A collisions, except 
for the data in Pb-Pb collision, where the hadronic absorption mechanisms failed to follow 
the abnormally strong J/ip suppression. 

The outline of the paper is as follows: In section II, we describe the transport approach 
and cascade method. In section III, the nucleus-nucleus interaction in GCP is described. 
In section IV, the nucleon and comover absorption mechanisms for J ftp suppression is 
investigated. In the last section, some further discussions are given. 

II. Transport approach and cascade method 

In a many body system, if correlation between particles is weak, i.e. independent 
particle pictures is a good approximation, and if the interaction is dominated by short- 
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range forces, we may treat the system as a collection of point particles propagating freely 
between successive short-range interactions. This system can be described by a set of 
Boltzmann equations, which can be solved using the cascade algorithm |]1T| |T^. 



In the absence of interactions, particles travel on straight line trajectories, and the 
equation for the time evolution of Ua is 

p^df,UJa{x,t,P) =0. (1) 

uja is the probability density of finding an on-shell particle a, with momentum P, at 
position = {x,t). Considering the interaction between particles, an S"— matrix is 
included to determine the interaction going from the initial state \bi,b2, ■ ■ ■ ,bn > to the 
final state |ci, C2, ■ ■ • , >, 

n m 

< ci, C2, ■ ■ ■ ,Cm\S\bi, 62, ■ ■ ■ >= A„^„(27r)^5^(E Pfc, - EPcJ- (2) 

i=l j=l 

Summing over all possible incoming channels and outgoing channels, we get a very general 
transport equation, with local interactions, in the absence of mean fields, 

d^Pt 



p^d^,uJaix,t,V) =Y. E j n {2^Y2Eb ^''^(^'^'^ 



b\,b2,--,b 



E E J n (2Tr)32E 

n m 

(2vr)V(Ep6,-EPcJ 

1=1 k=l 

n m 

[- E ^ab.6\p - fb.)2Eb, + E Sac,5Hp - fc,)2E,,l (3) 
i=l j=l 

where '— ' and '+' in the fourth line before 5 functions represent the decrease of particles 
of type a from the incoming channels and the increase of a in the outgoing channels, 
respectively. In principle, given the initial values of the distribution and all A^^m's, 
Equation (Q) can be solved. 

A particle cascade is one of the best method for solving these equations. In a cascade, 
the probability distributions Ua of particle a are sampled by Ma test points, 

^ ^ N ^ '^^ ^ ^ ^ 

COa{x,t,p) = J^E^'i^ - ^^ ii))S'iP - P^)' (4) 
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where Na is the total number of a. In the hmit Ma ^ oo, uo a can be described accurately 
by the density of test points. In the collision free limit, 

Ma j=i 

is a solution to Eq.(|l]), where tJi=Pi / Ei. Between collisions particles travel along a 
straight line, as in Eq.(|). 

For the case of resonance decay (n=l), 

m 
k=l 

= Vb^Cl+C2 + - + Cm- (6) 

?7b_»ci+c2+ --+c™ represents the decay rate for the channel of b Ci + C2 + ■ ■ ■ + Cm- The 
total decay rate is rj^ = J2 J2 Vb^ci+c2+---+cm- can convert the decay rate to the 
resonance lifetime, 77^ = l/(7brfe), where 7^ = \ Ph IZ-E'b is the Lorentz decay factor, is 
the inherent lifetime of resonance state h. 
In two body collisions (n = 2), 



1 (fp^ 



■ 2Et2 J /ii (27r)32E, 



2— >m 



|2 



■3 

m 



(2vr)V(p,,+p,,-EPcJ 

k=l 

= \V2- Vi |(T5^+52^ci+C2 + -+c™, (7) 

where ai,^+i,^^ci+c2+---+cm is the partial cross-section for the channel of 61 + 62 — ^ Ci + C2 + 
■ ■ ■ + Cm- The total cross-section is found by summing over all outgoing channels, 

^b^+b2 =51 X! '^6i+62^ci+C2+--+c„- (8) 

m Cl,C2,---,Cm 

The product | ^2 — "^i | • crbi+b2 is a cylindrical volume with a length | — -wi | and a 
cross-section ah^+b2- K represents the probability of collision between particles hi and 62 
in a unit time. We can simulate this term by making a collision between a point particle 
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bi and a point particle 62 whenever they are approaching each other within a cross-section 
cr^i+bj. We can set the coUision time to be the time when the distance between the two 
particles is at a minimum. When a collision occurs, the partial cross-sections are used to 
determine the branching ratio to a particular channel. 

The momentum distribution for the outgoing particles are selected within the phase 
space, weighted by and constrained by the energy momentum conservation ||16[| , 

TT ^ \A |2 

(2vr)32E,J^"' 

n m 
1=1 k=l 

(9) 

The cascade reduces the solution of the transport equations into the scattering of a set of 
classical point particles with known cross-sections and branching ratios. However, there 
is one problem unsolved. Although the transport equation is local and Lorentz invariant, 
the cascade breaks this invariance by allowing particles to collide at a distance d = y^cr/vr 
apart (for n = 2). The collisions are time ordered and this ordering is frame dependent. 
Therefore, Lorentz invariant is not strictly satisfied. The codes based on cascade method 
generally have this problem. The result of Monte Carlo simulation depends on the se- 
lection of frame more or less. At present energy range, the results of different frames 
consistent with each other qualitatively. However, the problem will be more serious at 
RHIC energy region. In ref. |TT| a method is proposed to recover this invariance. 



III. The nucleus-nucleus interaction in GCP 

In building the interaction of GCP, it is convenient and instructive to build up the 
physics step by step. 

1. Nucleon-nucleon interaction 

In GCP, the processes of the resonance production, vr production and the resonance 
decay are included in nucleon-nucleon interaction table, 

N + N — ^N* + N*, (10) 
A^ + A^ — >N + N + lTr, (11) 
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N* — ,N + 7i, (12) 

where I is the number of produced tt, which is related to the energy in the center of mass 
frame. The probabihty of the interactions is determined by the branching ratios. In GCP, 
the difference between proton(p) and neutron(n) is not considered. 7r+,7r^,7r° are also 
regarded the same as tt. 

(1). The process of vr production 

According to the experimental results, the average number of produced charged par- 



ticles in N-N collisions can be expressed as ||T7| 

< Nch >= 0.88 + 0.44/ns + 0.118(/ns)^ (13) 

where s is the square of the energy in the center of mass frame, in unit GeV^. Since 
the difference of tt^,tt^ and tt" is neglected, after removing the contribution of leading 
particles, we can get the average multiphcity of pions 

<n>=[< > -1.5] * I (14) 

Introducing the selected multiplicity distribution, one can use GCP to get the multiplicity 
n for each collision. The energy and momentum conservation are constrained during the 
collision processes. No other dynamics mechanism is included in GCP. In GCP, we choose 
Pj-limited phase space distribution. Considering the leading particle effect of incoming 
nucleon, uniform longitudinal momentum distribution of pions is selected. We put KNO 
distribution into GCP and get the multiplicity distribution of pions which is shown in 
fig.l. The statistical results of the rapidity distribution and the transverse momentum 
distribution of the produced particles in N-N collisions is shown in fig. 2. 
(2). The process of resonance decay 

A^* is a resonance state which is produced during N-N interaction. In GCP, A^* will 
decay to N and vr within a given time r, which is shown in Eq.(0). 
2. p-A and A- A collision 

The description that vr is produced immediately after each N-N collision in p-A and 
A-A collisions is not actually the case. The longitudinally excited wounded nucleon 
model(LEWNM) is introduced in GCP. In LEWNM, each N-N colhsion leads to two 
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longitudinally excited clusters. The clusters decay after a certain time. The decay of 
clusters is the source of other produced particles in N-N coUision. 

In LEWNM the nucleons are not transversely excited but only longitudinally stretched 
after the N-N collision. The collision of a beam nucleon and a target nucleon will result 
in two excited clusters, 'string!' and 'string2': 

N + N — > stringi + string2- (15) 

The mass of the clusters are fixed in the following way: By using GCP, the process of 
TT production in N-N collisions is simulated. In the center of mass frame, the particles 
are divided into two clusters according to their z direction. The invariant masses of all 
particles in each of the two clusters are counted over a broad range of cm. energies of 
the N-N system. Then we can get the empirical formula of the mass of the cluster which 
is a function of the cm. energy of the N-N system. 

The longitudinally excited clusters 'stringi' and 'string2' subsequently decay, 

string — > N + It:. (16) 

If the resonance or string collide with other particles before decaying new excited clusters 
may produce. Decay of the clusters are the source of particle production. 

The secondary collisions are also considered in p-A and A-A collision. The processes 
mainly include: 

N {N*, string) + 7r — >N + Ti + l7:, (17) 

TT + TT )• TT -|- TT -|- ZtT. (18) 

The process of N-N collision is not changed after introducing the LEWNM into GCP. Ac- 
cording to the experimental condition, wc also have simulated the S-U collision. The sta- 
tistical results of transverse energy(£'r) distribution are shown in figure 3. The LEWNM 
gives a good description of many observable quantities in N-N, p-A and A-A collisions. 
It can be applied to the study of relativistic nucleus-nucleus reactions. 
3. CoUision geometry 
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The nuclear geometry plays an important role in relativistic heavy ion collisions. The 
overlapping area of the two colliding nuclei is determined by the impact parameter b. We 
call the nucleons in the overlapping area the participants. The number of participants 
is calculated according to the nuclear geometry. Considering the change of b, we get 
the distribution of the number of participants. The larger the b, the less the number of 
participants. In GCP, b is an input parameter. It can be fixed or selected randomly in a 
certain region. 

In section II, we have discussed that the collision is possible when the distance of 
closest approach of the two colliding particles is less than the interaction range d, related 
to the total cross-section hy d — ^Jojn. The possible future coUisions are ordered in time 
and form the collision list. The next collision is the earliest one on the list. 

So the basic considerations of GCP only relate to coUision geometry and cascade 
simulation, no other dynamics of the production process is included. The aspects we 
mentioned above are only related to the nuclear geometry and kinematics. The detailed 
physical assumptions can be introduced according to the physical models. This is the 
remarkable advantage of GCP. 

IV. The hardronic mechanisms on J/ if; suppression 

Two hardronic mechanisms of J/ip suppression are added to GCP to simulate the 
production and absorption of J/ip. 

1. The process of J/ip production and absorption 

The channels for J /if; production and absorption are added to GCP to simulate the 
process of J/%1) production and absorption: 

N{N\ string) + N{N\ string) — > N + N + J/^ (19) 
N{N*, string) + J/ijj — >N + D + D (20) 
TT + J/^ — ^ D + D (21) 
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Since the probability of J/ip production in N-N collision is relatively small, the cross- 
section of J / ip production is enlarged in the simulation to increase J / ip production prob- 
ability. In general, one or two J/ip are produced in each event in the simulation. 

2. The mechanism of the J/ip absorption in nuclear matter 

Based on above discussions, one can accept the following physical picture that the 
J/ip production can be divided into two steps. The first step is the production of a cc 
pair, which is produced perturbatively and almost instantaneously. The second step is the 
formation of a physical state of J /ip., that needs a much longer time. In a nucleus-nucleus 
collision, cc are produced by hard scattering processes of beam nucleon and one target 
nucleon. The produced cc may interact with another target nucleon and these cc-nucleon 
interactions may lead to the break up of cc via the reaction 

CC + N{N*, string) — > D + D + x, (22) 

which turns cc into DD pair. Thus, cc-nucleon interactions will give a suppression of J/ip 
production. 

3. The mechanism of the interaction between J/ip and produced mesons 
Comovers usually refer to the secondaries produced in high energy heavy ion collisions, 

such as 71, p and u mesons, etc. In A-A collisions, besides J/ T/^-nucleon absorption, J/ ip 
particles also suffer interaction with secondaries that happen to travel along with them, 
which also causes J/ip suppression. The reaction is shown in Eq.([2T|) . 

The details of simulating the process of J / ip production and absorption are as follow- 
ing: 

(1) . An entry for J/ip production is added to the interaction table. We use average 
impact parameter < 6 > to discuss minimum biased data. For the situation of different 
Et bin the values of b are given by experimental groups f^. In order to increase the 
probability of Jjip production, we enlarge the cross-section of J/ip production to ensure 
one or two J /ip produced in every event. Running GCP, we get statistical number of J /ip, 
Npro, in 10000 events. 

(2) . Put the channel for J/?/'- nucleon interaction into GCP. Using the same condition 
as (1), we run GCP again to count the number of J/ip, Nabs{N), in 10000 events, where 
the nucleon absorption is included. 
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(3). At last, J / ip-comover interaction is added to GCP. Using the same condition as 
(1), we obtain the number of J/ip, Nabs{N+co), after considering the nucleon and comover 
absorption. 

The J/ip survival probability in A-B collision is expressed as 

S = Nabs/Npro. (23) 

The absorption cross sections of J/ip, i.e. the cross-sections that the produced J/ip 
particles interact with target nucleons or produced secondary particles, are important 
parameters to explain the J/ip suppression in hardronic environment. There are two pa- 
rameters, J/'?/'-nucleon cross-section aabs{N) and J/'i/'-comover cross-section <7abs{co)- After 
analyzing many sets of experimental data in p-A collisions, Gerchel and Hiifner fl^ found 



that the experimental J/ip production data for p-A collisions can be fitted well with an 
effective J/ip-imdeon cross section a = 6.2mb or 6.9mb. To account for the A-A data, 
J/V'-comover absorption cross-section is generally regarded as about 3mb. In GCP, cTabs{N) 
and (Tabs(co) are input parameters. We wish to adjust the parameters to agree well with 
the experimental data oi J /ip suppression. 

First we discuss the case of minimum biased data. Considering only the absorption 
of J/ip by nucleus, the GCP simulation of Jjil) survival probabilities at SPS energy are 
expressed as squares in fig. 4. The experimental data are shown as black triangles. The 
J/'i/'-nucleon cross-section is taken to be aabs{N) = 7m6, which is in good agreement to 
earlier works |T^. The lifetime of the wounded nucleon is taken to be Ifm/c. One can 



see clearly that our simulation is in good agreement with the experiment data of J/ip 
production in p-A collisions. Most of the J/ ip suppression data for A-A collisions are also 
fitted, but one can not explain the data in Pb-Pb collision. Using GCP, we can simulate 
the J/ ip production and absorption in p-A and A-A collisions up to S-U data successfully. 
Next, J / ip-comover interaction is added to GCP. Considering the nucleon and comover 
absorptions of J/ip, we repeat the above procedure again. The value of crabs{N) is taken the 
same as above. The J / ip-comovei absorption cross-section is <7abs{co) = 2m6. The results 
are expressed as open triangles in fig. 4. It shows that the J/ip-comovei interaction in A-A 
collisions is more important than that in p-A collisions. However, to include the comover 
absorption has not introduced qualitative difference from the result including only nucleon 
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absorption. The combination of these two absorption mechanisms of J/%1) still can not 
explain the data for Pb-Pb coUision. In the above situation, the lifetime of the wounded 
nucleon is taken as Ifm/c. In GCP, the decay of the wounded nucleons is the source 
of particle production. If their lifetime decreases, other particles such as mesons will be 
produced earlier. This may increase the probability of J/-(/^-comover interaction. Fig. 5 
shows the results in p-A and A-A collisions obtained using the lifetime of the wounded 
nucleons as 0.2fm/c. It can be seen from the result that the shorter lifetime does provide 
a much stronger comover absorption. Now the Pb-Pb data can be reached. However the 
other fitting points are much lower than the experimental data in lighter A-A collisions. 
The above analysis tell us that one can not explain the J/ ^0 suppression for all the observed 
experimental data consistently based on the one set of parameters. Then we adjust crabs{N) 
to 6mb, aabs{co) to 3mb, repeat the above step again. In this case, the comover absorption 
of J/ip is more manifest. But the overall results are similar to fig.4. 

Now we turn to discuss the result of different Et bins. In the same way as above, we 
simulate the two kinds of J/ip suppression, the interaction of J/ip with nuclear matter 
and produced meson in S-U and Pb-Pb collisions. For different Et bins, the values of 
b are given by experimental groups, which are listed in table 1. The parameters aabs{N) 
and aabs(co) used in considering the J/'^-nucleon absorption and J / ip-como\er absorption 
and the hfetime of the wounded nucleon are taken to be the same as those used for fig.4. 
in fitting the minimum biased data. The squares in fig. 6 show that J/ip suppression in 
S-U collision could be described very well based on our nucleon absorption. However, 
the fitting for the last E'r-bin in Pb-Pb collision using only nucleon absorption shows 
that (the most right square in fig. 6) the data in Pb-Pb collision could not be fitted using 
the same set of parameters. The interaction of J/i(j with produced mesons make further 
suppression of J/ip in S-U and Pb-Pb collisions, which are shown as triangles. But the 
combination of these two kinds of absorption (see "y " in fig. 6) still can not explain the 
extra strong J/ip suppression in Pb-Pb collision. There may be anomalous suppression 
caused by other mechanisms, e.g. QGP formation. 



V. Results and Discussions 
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GCP is designed such that the cascade algorithm is separated from the physics models. 
In this way the problem of the cascade can be isolated and studied separately, and different 
physics models can be compared in a common framework. GCP itself is only a tool to 
simulate the cascade process. The various physics models can be built based on it. This 
is the key difference between GCP and most other cascade codes. 

In this paper, we focus our attention on the Monte Carlo simulation of J /ijj produc- 
tion and absorption directly by using GCP. Our work is divided into two parts. First, the 
longitudinally excited wounded nucleon model(LEWNM) is added to GCP. It is used to 
handle the nucleus-nucleus interaction. The LEWNM is found to give a good description 
of many observables in N-N and p-A collisions, such as the multiplicity distributions, the 
rapidity distributions, and the transverse momentum distributions of the produced par- 
ticles. According to the experimental conditions, we have also simulated the transverse 
energy distribution for S-U collision. Next, two hardronic mechanisms of J/V' suppression, 
J/-?/'-nucleon interaction and J / ip-comoYei interaction, are added into GCP to simulate 
the production and absorption of J/ip in p-A and A- A collisions at SPS energy. First 
we consider the case of minimal biased data. Then the results of different Et bins are 
discussed. Using the nucleon absorption mechanism, the results agree well with the ex- 
perimental data of J/'^ production, except for those in Pb-Pb coUision. Based on above 
discussion, the comover absorption mechanism is added to GCP to study the suppression 
of J I'll) production. The J/tl) suppression in Pb-Pb collision can not be explained even 
by the combination of the above two mechanisms. It seems that some new mechanisms 
is needed to study the anomalous J/^ suppression in Pb-Pb collision. This may indicate 
the formation of QGP. 

Acknowledgments 

This work is supported in part by the National Natural Science Foundation of China 
and Natural Science Foundation of HeBei province. We would hke to thank Zhang X.F. 
and Tai A. for stimulating and helpful discussions. 



13 



REFERENCES 



[1] T. Matsui and H. Satz, Phys.Lett. B178 (1986) 416; T. Matsui, Z.Phys. 
038(1988)245. 

[2] The NA38 Collaboration, Z.Phys. 038(1988)17; Phys.Lett.B220(1988)471; 
B255(1991)255; Nucl.Phys. A544(1992)209c. 

[3] C. Morel et al., Phys. Lett. B 252 (1990) 505; D. M. Aldel et al. Phys. Rev. Lett. 
66 (1991) 133. 

[4] M. Gonin et al. (NA50), Report at Quark Matter 1996, Heidelberg, Germany; 

[5] F. Fleuret (NA50 Collab.), in Proceedings of the 32nd Rencontres de Moriond, QCD 
and High Energy Hadronic Interactions, Les Arcs, France, 1997. 

[6] D. Kharzeev, invited talk at QM'97. 

[7] J.-P.Blaizot and J.-Y. Ollitrauh, Phys. Rev. Lett. 77 (1996) 1703. 
S. Gavin and R. Vogt, Phys. Rev. Lett. 78(1997) 1006. 



[9] A. Tai, W. Q. Chao and X. X. Yao, |hep-ph/9701207 . 



[10] D. Kharzeev and H. Satz, Phys. Lett. B366 (1996) 316. D. Kharzeev, C. Lourenco, 
M.Nardi and H. Satz, |hep-ph/96122r7| . 



[11] Y. Pang, in Proceedings of CCAST Symposium/ Workshop on Particle Physics at 
Fermi Scale, Beijing, May 1993, Edited by Y.Pang, J.Qiu and Z.Qiu( Gordon Breach, 
New York, 1993), p451. Ion Collisions. 

[12] Gassing W, Bratkovskaya EL. Nucl.Phys.A623(1997)570; N.Armesto, A.Gapella, 
|hep-ph/ 97052751 . 



[13] B.H.Sai, A. Tai, |nucl-th/9803033 



[14] Y. Pang, RHIG Summer Study '96: Theory Workshop on Relativistic Heavy Ion 
Collisions July 8-19, 1996 at Brookhaven National Laboratory, edited by D.E.Kahana 
and Y.Pang, pl93. 



14 



[15] C. Gerchel and J. Hiifner, Phys. Lett. B207, 253(1988); C. Gerchel and J. Hiifner, 
Nucl. Phys. A544, 513c(1992). 

[16] Y. Pang, Lorentz Invariant Multiparticle Phase Space Besed on ARC Experience. 

[17] Cheuk-Yin Wong, Introduction to High Energy Heavy Ion CoUisions, world scientific, 
printed in Singapore. 



15 



Figure Caption 



Fig.l: the statistical results of the multiplicity distribution of pions in N-N coUision. 

Fig. 2: the statistical results of the rapidity distribution and the transverse momentum 
distribution of the produced particles in N-N collision. 

Fig. 3: the statistical results of the transverse energy distribution of the produced 
particles in S-U coUision. 

Fig.4: The J/-^ survival probability in p-A and A-A collisions obtained using the 
lifetime of the wounded nucleon as Ifm/c, together with minimum biased experimental 
data for(from left to right) p-Cu, p-W, p-U, 0-Cu, 0-U, S-U and Pb-Pb coUisions . 

Fig. 5: The same as Fig.4 using the hfetime of wounded nucleon as 0.2fm/c. 

Fig. 6: The J/'^ survival probabihty for different Et bins together with experimental 
data in S-U and Pb-Pb collisions. 



Table Caption 

Table I. The value of < 6 > for different Et bins in A-A collisions. 
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Table I 
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